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Abstract 

We present numerical results for non-compact three-dimensional QED for numbers of flavors 
Nf = 1 and Nf = 4. In particular, we address the issue of whether chiral symmetry is 
spontaneously broken in the continuum limit, and obtain a positive answer for Nf = 1, with 
a dimensionless condensate estimated to be /3 2 (ipif}) ~ O(10 -3 ), implying that the critical 
number of flavors Nf c > 1. We also compare the Nf = 1 and Nf = 4 models by analysing 
the transition from strong to weak coupling behaviour using an equation of state based on a 
continuous phase transition. While some qualitative differences emerge, it appears difficult 
to determine whether Nf = 4 lies above or below Nf c . 
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1 Introduction and Motivation 



The nature of the ground state of Quantum Electrodynamics in 2+1 dimensions and its 
dependence on the number of fermion species Nf has been an important problem for non- 
perturbative field theorists for 20 years pp. The question is the value of critical number of 
flavors Nf c such that for Nf < Nf c the theory's chiral symmetry is spontaneously broken, 
resulting in a dynamically generated mass for the fermion. Over the years the main method 
of attack has been to solve the Schwinger-Dyson equation for the fermion self-energy 
to see if limfc^o 7^ 0; this has resulted in a variety of predictions for Nf c ranging between 
32/7T 2 and 00, with most recent attempts yielding Nf c ~ 4 [2J. Over the same period there 
have also been attempts to decide the issue by numerical simulation of non-compact lattice 
QED [31I31E1- Numerical approaches are hindered by large finite volume effects, and the 
intrinsic smallness of the relevant signal, the so-called dimensionless condensate /3 2 (?/>?/>) (to 
be defined in what follows) in the continuum limit: in a recent study jB] of the model with 
Nf = 2 we placed an upper bound on this quantity of O(10~ 4 ), but even this does not 
exclude the predictions of some self-consistent approaches with Nf c > 2 [ZUHj- 

Interest in the problem has recently been revived by the suggestion that QED3 may be 
an effective theory for the underdoped and non-superconducting region of the phase diagram 
of high-T c superconducting cuprate compounds EH . In brief, superconductivity in these 
substances is confined to planes defined by Cu02 layers, thus motivating a 2+ld description. 
The superconducting order parameter has a <i-wave symmetry, implying that there are four 
nodes in the gap function A(k) as the Fermi surface (which in 2+ld is a curve) is circum- 
navigated. At each node the low-energy quasiparticle excitations obey an approximately 
linear dispersion relation E(k) oc \k — k no d e \, with the result that it is possible to rewrite the 
action for eight distinct low energy species (spin up and spin down at each of four nodes) in 
a relativistically invariant form in terms of Nf = 2 species of four-component Dirac spinor 1 . 
The next link in the chain is the observation that disruption of superconducting order in 
a planar system occurs via condensation of vortex singularities in the phase <fi of the Un- 
valued order parameter, while |A| remains unchanged. Because the order parameter field 
is doubly charged since it arises from electron pairing, in the presence of a distribution of 
vortices it is impossible in general to reabsorb (f) into the definition of the quasiparticle fields 
while leaving their wavefunction single-valued. There is, however, a particular gauge ^T] 
in which the problem can be recast in terms of single-valued quasiparticle fields interacting 
with a real-valued vector field. The component of the vector field that is minimally coupled 
to the quasiparticles remains massless under quantum corrections (vacuum polarisation due 
to virtual particle hole pairs), and is identified with the photon of QED3 2 . It is then ar- 
gued [5J HO] that the effective action for the photon fluctuations is of the non-compact form 
F^ u , with an effective "electric charge" coupling photons to quasiparticles proportional to 
the vortex disorder parameter, which is non-zero as soon as vortices unbind, ie. outside the 
superconducting region of the phase diagram. 

1 For phenomenologically relevant models the action for each individual flavor exhibits a spatial anisotropy, 
a feature ignored in this paper. 

2 In fact, the photon can be regarded as the Goldstone boson associated with vortex condensation |12|. 
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If QED 3 is a relevant effective theory for cuprates, then the abstract theoretical problem 
of the value of Nf c assumes concrete phenomenological reference. If Nj c > 2, then the 
theory is chirally broken at zero temperature. On retranslating from the Dirac spinor basis 
to the original electron degrees of freedom, the chiral order parameter is reinterpreted as 
an order parameter for spin density waves, whose wavevector gets shorter and shorter as 
doping is decreased, until at zero doping the Neel antiferromagnetic state is recovered |lUj . 
This picture therefore predicts the existence of a phase boundary between superconducting 
(dSC) and antiferromagnetic (AFM) phases at some non-zero doping in the zero temperature 
limit. If, on the other hand, Nf c < 2, the chirally symmetric ground state manifests itself 
as a tongue of "pseudogap" phase separating dSC from AFM, in which normal Fermi liquid 
properties may be modified as a result of a non-perturbative anomalous dimension for the 
fermion field ||. 

In order to motivate further study, let us review two arguments concerning Nf c which 
are independent of Schwinger- Dyson analyses. A recent discussion of similar issues can be 
found in [7j. The first due, to Appelquist, Cohen and Schmaltz [T3], postulates an inequality 

flR < fuv (1) 

between thermodynamic free energy densities (ie. pressures) due to the weakly-interacting 
degrees of freedom characteristic of respectively long and short distance regimes. Since / 
is simply related to the number of weakly interacting degrees of freedom, this can be easily 
checked for QED 3 . In the UV limit fuv oc 1 + |(4iVy), where 1 counts the photon, which 
has only one physical polarisation in 2+ Id, and the second term, which includes a numerical 
factor arising from Fermi-Dirac statistics, counts free fermions. In the IR limit for Nf < Nf c 
there is still a photon, but chiral symmetry breaking U(2iV/) — > U(iV/) ® U(iV/) implies 2N 2 
Goldstone modes, so fm oc 1 + 2Nj. Eqn. (JTJ) then yields the bound Nf c < |. It should 
be pointed out, however, that this argument may not be applicable to a 2+ lei system, since 
in this low dimensionality bosonic modes exhibit strong fluctuations and hence are never 
weakly interacting |14j . 

Another argument starts from the similarity between QED3 and the 2+1 dimensional 
Thirring model defined by the Lagrangian 

G 2 

C-Thir = + m)tpi + ^-r{^il^i) 2 - (2) 

In the limit G 2 — > 00 there is a massless vector meson, ie. a ipip bound state with the same 
quantum numbers as the photon; moreover both models have identical quantum corrections 
when calculated in an expansion in powers of Ny 1 jT^j. However, the resemblance may also 
be non-perturbative in Ny 1 ; Itoh et al pointed out that in a certain non-local gauge the 
self-consistent gap equation for chiral symmetry breaking in the Thirring model coincides 
in the limit G 2 — > 00 with that of QED3, implying a vanishing chiral order parameter for 
some critical Nf c . In the Thirring case, however, the vanishing of the gap also happens for 
Nf < N'f c , G 2 = G 2 < 00, and signals a UV fixed point of the renormalisation group at 
which a non-trivial continuum limit may be taken. This picture has received support from 
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numerical simulations with finite Nf, which suggest that 4 < Nf c < 6 fTf\ ITS} IT5]. Since the 
pattern of global symmetry breaking is the same for both models, it is tempting to postulate 
that the IR fixed point behaviour of QED 3 coincides with the UV behaviour of Thirring in 
the G 2 — > oo limit, and hence Nf c = Nf c ~ 5. 

Since neither argument is obviously nonsense, and since they yield predictions for Nf c 
which differ significantly both from each other and from the Schwinger-Dyson approach, it 
is clear that consensus on this issue has yet to emerge, and that further numerical efforts 
are justified. In jB] we were unable to find any signal for chiral symmetry breaking, but 
also unfortunately were not able to definitively exclude it being present with the tiny value 
predicted in the self-consistent approach |H|. In this paper we build on that work by 
performing extensive simulations of non-compact QED 3 with Nf = 1 and Nf = 4, the hope 
being that if these two values were to straddle Nf c , then some difference in their behaviour, 
either quantitative or qualitative, may show up. The lattice formulation and simulation 
method is reviewed in the next section, and the numerical results presented in Sec. EJ We 
shall see that for Nf = 1 there is convincing evidence for broken chiral symmetry, implying 
Nf c > 1; unfortunately, while we can uncover some qualitative differences between Nf = 1 
and Nf = 4, there is no decisive evidence that they have different ground states in the 
continuum limit. 

2 The Model and the Simulation 

The action of the lattice model we study is 

8 N 

s = ^e p (i)e F (i) + ^^iWM(^')xiM (3) 

Z x,ii<u i=l x ,x' 

^/IV\J') = @X[l ^X~\-fl,l> ^X+U,fJL @XV 

M(x,x') = m 6 x ^ + -^rip(x)[6j^pJJ X p-6tfp-fiUl_fr lt ]. 

This describes interactions between N flavors of Grassmann-valued staggered fermion fields 
X, X defined on the sites a; of a three-dimensional cubic lattice, and real photon fields 9 XI1 
defined on the link between nearest neighbour sites x, x + jl. Since 2 is unbounded from 
above, Q defines a non-compact formulation of QED; note however that to ensure lo- 
cal gauge invariance the fermion-photon interaction is encoded via the compact connection 
Ux^i = exp(i8 x ^), with U x +ji-^ = U*. In the fermion kinetic matrix M the Kawamoto-Smit 
phases 

^(x) = (-1)^+-+^-! (4) 

are designed to ensure relativistic covariance in the continuum limit, and mo is the bare 
fermion mass. 

If the physical lattice spacing is denoted a, then in the continuum limit ad — > 0, (0) can 
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be shown to be equivalent up to terms of 0(a 2 ) to 



S = Y1 ^ j Md, + igAJ + mW + jF, v F, v (5) 

ie. to continuum QED in 2+1 euclidean dimensions, with ip, ip describing Nf flavors of 4 
component Dirac spinor acted on by 4x4 matrices 7^, and Nf = 2N. The continuum photon 
field is related to the lattice field via 8 xtl = agA^x), with dimensional coupling strength g 
given by g 2 = (a/?) -1 , and the field strength F M „ = d^Ay — d^A u . The continuum limit is 
thus taken when the dimensionless inverse coupling (3 — > 00. 

As reviewed in jB], for a > in the chiral limit the lattice action Q retains only a 
remnant of the U(2iVy) global symmetry of (0) under global chiral/flavor rotations, namely 
a U(iV) <g> U(iV) symmetry which is broken to U(iV) either explicitly by the bare mass 
m ^ 0, or spontaneously by a chiral condensate (xx) 0, in which case the spectrum 
contains N 2 exact Goldstone modes. It is expected that the symmetry breaking pattern 
XJ(2Nf) — > U(Nf) ® U(iVf) is restored in the continuum limit, implying the existence of an 
additional 7N 2 approximate Goldstone modes whose masses vanish as (3 — > 00. 

The simulations in this paper were performed using two different algorithms. The model 
with Nf = 4 is studied by exact simulations of the action with N = 2 using the hybrid 
Monte Carlo algorithm, exactly as done for Nf = 2 in [6 . The model with Nf = 1 is simu- 
lated using det^M as a functional measure using a hybrid molecular dynamics R algorithm 
[2*U] , which generates a Markov chain of representative configurations by evolution of a sys- 
tem of stochastic differential equations with timestep At. The algorithm has a systematic 
error oc N 2 At 2 ; in our work we use values of At ranging from 0.1 down to 0.0025, the value 
required getting smaller as (3 decreases, L s increases, and m decreases, and have checked 
that this is sufficient to render systematic errors smaller than statistical ones. Beyond this 
technical point, however, we should point out the conceptual issue that simulations with 
fractional N do not reproduce the physics of a local fermion bilinear action except in the 
deep continuum limit when flavor symmetries are restored. It remains an open issue within 
the lattice QCD community whether this is a source of significant systematic error at achiev- 
able lattice spacings [21]. Whilst we make no attempt to address this issue in the current 
study, we should not exclude the possibility that delayed flavor symmetry restoration could 
lead to a misleading prediction for the ground state properties at finite /3 [§]. In Sec. 13.21 
below we shall present evidence for flavor symmetry violating effects of just O(10%) at the 
weakest couplings examined. 

To check if chiral symmetry is broken in continuum QED 3 , we need to establish whether, 
on a system of finite extent L = L s x a 

Jim lim lim (3 2 (xx) ^ 0. (6) 

The order of limits is important: since a continuous symmetry never breaks spontaneously 
on a finite volume we must always work with rn^O, and only attempt to take the chiral limit 
once the thermodynamic limit is reached. Since for an asymptotically-free theory like QED3 
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the UV behaviour is governed by the gaussian fixed point at the orgin, all physical quantities 
are expressible in terms of the scale set by g, and hence to compare data taken at different 
lattice spacings, it is natural to use the combinations /3m , L s / (3 and (3 2 (xx) corresponding 
to the dimensionless combinations m/g 2 , Lg 2 , and g~ A (ipip) . As the continuum limit is 
approached data taken at different (3 should collapse onto a universal curve when plotted in 
these units. 

Taking all three limits of equation (jUJ) is extremely demanding computationally, as re- 
vealed in , where in practise we were only able to place an upper bound on the condensate 
for Nf = 2 of (3 2 (xx) < 5 x 10 -5 . The small size of the dimensionless condensate in this 
case implies that a large separation of scales £ ^> a may be taking place, and therefore a 
much larger infra-red cutoff L ^> £ required to probe chiral symmetry breaking if it is indeed 
present [22J. The small number is not, however, inconsistent with current estimates based 
on solutions of the Schwinger- Dyson equations |S] . It therefore appears unlikely that the 
critical value Nf c can be established using current computational resources if Nf c > 2. 

In this paper our aims are twofold. First, since the self-consistent solutions suggest an 
exponential suppression (ipip) ~ exp(— f(Nf)), where / is an increasing function of Nf for 
Nf < Nf c , the dimensionless condensate for Nf = 1 may well be much larger and hence 
easier to measure using current resources. Secondly, we will look for qualitative differences 
between data taken with Nf = 1 and Nf = 4 to see what evidence if any can be found 
that they lie on opposite sides of Nf c , by fitting the condensate data to a global equation of 
state based on a chiral symmetry restoring transition at finite (3 c (Nf). For Nf < Nf c this 
"critical" coupling need not correspond to a true phase transition, but instead could mark a 
crossover from strong to weak-coupling behaviour. 

3 Simulation Results 

For N f = 1 we have results from lattice volumes 24 3 (0.25 < (3 < 0.45), 32 3 (0.45 < (3 < 0.65), 
44 3 {p = 0.55), 48 3 {(3 = 0.65), and in addition five studies at (3 = 0.90 on 16 3 , 24 3 , 36 3 , 
54 3 and 80 3 . For N f = 4 we have 16 3 (0.15 < (3 < 0.2), 24 3 (0.15 < (3 < 0.225), 32 3 
(0.4 <{3< 0.6), and 48 3 ((3 = 0.4). For the most part we examined mass values in the range 
0.003 < uiq < 0.025. For each parameter set we took typically 200 - 300 trajectories of mean 
length f = 1.0. Several different computing facilities were used, so an estimate of the cpu 
resource required is approximate at best, but a figure of 12 months running on a cluster of 
30 1.2 GHz processors is not unreasonable. In Figs. Q (Nf = 1) and|2] [Nf = 4) we show 
the datasets thus obtained for (xx) as a function of f3 and mo- The lines are fits to a global 
equation of state to be described in Sec. 13.31 

3.1 Chiral Condensate in Chiral and Continuum Limits 

In an attempt to get an overview of the model's behaviour, for each (3 and Nf we first 
attempt an extrapolation to the chiral limit of the form 

3 

(XX) = (xx)o + Am + Bml (7) 
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and plot the resulting (xx)o as a function of (3 in Fig. El We have checked that at least 
at stronger couplings where the signal is appreciable the results are relatively insensitive to 
changes in the extrapolation function, eg. to incorporate non-analyticities due to Goldstone 
degrees of freedom We see that for both Nf values chiral symmetry clearly appears 

to be broken at strong coupling, but that (xx)o decreases as j3 is increased, so that we can 
identify either a crossover or a symmetry restoring phase transition at (3 ~ 0.5 (Nf = 1) or 
(3 ~ 0.25 (Nf = 4). Global fits to the data will be used in Sec. 13.31 to explore the nature 
of these points further. There is some evidence already from Fig. |3] that the curvature of 
(xx)o(P) differs between the two values of Nf, a feature already noted in j^j. 

In an attempt to see whether chiral symmetry is broken in the continuum limit (3 — > oo 
we compare plots of the dimensionless variables (3 2 {xx) versus (3rriQ for different Nf in Fig. 0] 
The data for Nf = 2 are from [(3,, where it was concluded that when extrapolated to the 
continuum limit /3 2 (xx) ~ 10 -4 . By comparing Nf — 1, 2 and 4 one sees that increasing Nf 
at fixed (3mo suppresses the condensate, as might be expected due to enhanced screening 
from virtual XX pairs. There is clearly no prospect of finding a signal for chiral symmetry 
breaking with Nf = 4 with this dataset. The Nf = 1 data, however, are significantly larger 
near the chiral limit. 

In Fig. |5]we show the plot for Nf = 1 zoomed in the neigbourhood of the origin for two 
different lattice spacings each with two volumes. The smaller of the two (3 = 0.90 lattices is 
intermediate in physical volume between the two (3 = 0.65 lattices, whereas the (3 = 0.90 80 3 
data is closest to both continuum and thermodynamic limits. There are clear finite volume 
effects at both lattice spacings for /3m < 0.005. By comparing (3 = 0.65 with f3 = 0.90 
it is clear there are also still lattice discretisation artifacts of O(10%). A linear fit to data 
on the larger volume extrapolates to a chiral limit value (3 2 {xx) = 0.0041(4) for (3 = 0.65 
(fit to the lowest 5 mass points), and 0.0043(2) for (3 = 0.90, with the fitted slopes differing 
by about 20% 3 . Due to finite volume effects these numbers are if anything underestimates. 
Therefore whilst it is premature to quote a value for the dimensionless condensate due to 
lack of control over systematic errors, we are confident in claiming that it is ^ 10~ 3 , ie. 
greater than zero, almost two orders of magnitude greater than the upper bound quoted 
for Nf = 2 in [S], and very roughly half the result for Nf = (ie. quenched QED3) found 
in j^]. Note, moreover, that our lowest mass datapoint at (3 = 0.90 has (3 2 {xx) oru y three 
times its value in the chiral limit, making the linear extrapolation of Fig. El considerably less 
adventurous than those of j3] . Any extrapolation which obtains a significantly smaller result 

in the chiral limit must be based on an analytic form which curves sharply as the chiral limit 

1 

is approached. One possible form is (xx) — a o + Oi"t.q H , predicted by the effective chiral 

theory describing the influence of pion fluctuations (see eg. [23 ). The applicability of such 
a model, however, presupposes chiral symmetry breaking. Our conclusion that Nf c > 1 is 
probably quite robust. 

3 Note that the slope of the line corresponds to the longitudinal susceptibility xi — d 2 InZ/dni^, which as 
a two-point function is expected to be more severely finite-volume affected. 
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3.2 Meson Spectroscopy 

A physically meaningful way of assessing how far the simulations lie from the continuum 
limit comes by checking the degree of flavor symmetry restoration in the meson spectrum. 
Recall that in the continuum theory with Nf = 2 [21], chiral symmetry breaking U(4) — ► 
U(2)(g)U(2) results in eight Goldstone bosons, which somewhat sloppily we will continue to 
refer to as "pions". If we denote the Euclidean time direction as x%, then the pion states can 
be divided into 4 pseudoscalars $(75 $(75 ®Ti)ijj and 4 scalars ^(74 ^(74<8>Tj)^, 

where {H.2,t} generate a U(2) symmetry acting on flavor indices. In the lattice model the 
corresponding states in terms of staggered fermions read 



^(75 & x(x)x(x)e(x); (8) 

^(74®r M )^ ^> + fi)e{x)C^x)] (9) 

^(75 ® t m )V> ^ I^aIxWx^ + i> + \)e(x)C v (x)Cx(x); (10) 

^(74®!)^ |^a|x(s)x(^ + £ + z> + A), (11) 

where we have defined additional phase factors e(x) = (— l^i+^a+xa^ ^(x) = e(a;)7/ At (a:)(— l)^, 
and an average over forward and backward link shifts is understood. The parity operation 
for staggered fermions with the definition (@J) for the Kawamoto-Smit phases is 

x = (x 1} x 2 ,x 3 )i-+a/ = (l-xi,x 2 ,x 3 ) ; x(x) ^ e(x')x(x') , x(x) ^ e(x')x(x f ) (12) 

The assignment of scalar and pseudoscalar labels to the states is readily checked. 

Spectroscopy on the lattice through decay of temporal correlation functions is most trans- 
parent when restricted to operators defined on a single timeslice. This means that in addition 
to the state (JBJ) which interpolates an exact Goldstone state for the residual lattice chiral 
symmetry breaking U(l)(g)U(l)— > U(l), it is possible to examine two scalar states interpo- 
lated by one-link operators (jHJ) and one pseudoscalar interpolated by a two-link operator 
(fTUJ) . At non-zero lattice spacing these latter states are not associated with Goldstones of 
an exact symmetry, and the degeneracy between the states expected in the continuum limit 
is split by terms which are generically 0(a). The extent to which the measured masses are 
degenerate, therefore, gives some measure of the approach to the continuum limit. Note that 
flavor symmetry restoration should still manifest itself in this fashion even if chiral symmetry 
is unbroken, since in this case the various states are related by U(2)®U(2) rotations. 

We have tested the approach to the continuum limit for Nf = 1 by measuring the masses 
of the local, 1-link and 2-link pions at (5 — 0.6, 0.75 and 0.9 with fixed values of f3rriQ = 0.0009 
and L s /j3 = 40. In this way any systematic errors due to eg. finite volume should cancel. 
The results are plotted in Fig. |U1 Whilst the data is still too noisy to determine the level 
ordering unequivocally, we can see that the splitting of 40% at (5 = 0.6 is reduced to 15% by 
(3 = 0.9, which thus gives a sample estimate for the residual systematic errors due to lattice 
artifacts at this closest approach to the continnuum limit. 
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3.3 Global Fits to Equation of State 

In this section we change tack, retreating from direct attempts to probe the continuum limit, 
and instead exploiting the data at stronger couplings to make qualitative statements about 
the behaviour as a function of Nf. In the strong coupling limit (3 — > 0, it is known rigorously 
that chiral symmetry is spontaneously broken and indeed our data confirm this (see 
Figs. I1I2|) . Therefore as (3 increases, for Nf > Nf c there must be a chiral symmetry restoring 
phase transition for some finite (3 C . For Nf < Nf c , since we believe the order parameter 
is exponentially small in the continuum limit, the relic of this transition may persist as a 
crossover between strong and weak coupling behaviour, much as in lattice QCD. For Nf 
exactly equal to Nf c , the order parameter should vanish in the chiral continuum limit, but 
with a non- analytic response to a small bare mass of the form (xx) °^ m oi a behaviour 
characteristic of a so-called "conformal fixed point" fI7\ I24j . 

We therefore proceed with the assumption (which must be unphysical for Nf < Nf c ) that 
there is a continuous chiral symmetry-restoring phase transition described by the following 
equation of state: 

m = A(P-p c )(xx) p + B(xx) 5 , (13) 

and attempt to fit the entire dataset, or at least a significant fraction, for each Nf. Eqn. 
has 5 free parameters, the "critical" coupling j3 c , exponents 6 and p, and the amplitudes A 
and B. Setting f3 = /3 C shows that 5 has its conventional interpretation as a critical exponent 
describing response of an order parameter to an external field, and setting mo = shows that 
the conventional exponent describing scaling of the order parameter in the broken symmetry 
phase is given by j3 mag = (5 — p)~ l . The fits are performed using the MINUIT package. 

For Nf = 1 there are a total of 151 points in the dataset (excluding (3 = 0.9). The best 
fit found was to the 109 points with 0.25 < (3 < 0.45 using data from all masses and from 
the largest available volume at any given (3. This gives 

A = 4.43(17) 5 = 2.33(4) p c = 0.451(4) X 2 /dof = 1.02 , , 

B = 1.87(7) v = 1-50(2) ^ > 

Acceptable fits persist, with little change in the parameters, if data with (3 = 0.400, 0.450 is 
excluded. The x 2 /dof increases to 1.5 if {3 = 0.55 is included, and 2.5 if (3 — 0.65 is included, 
once again with relatively little impact on the fit parameters. If instead strong coupling data 
with (3 < 0.40 are excluded, then x 2 /dof rises to ^ 4, and the fitted values of 5 and p drift 
downwards. 

The raw data are plotted with the fit (|14|) in Fig. Q and the fit is also shown in "Fisher 
coordinates" for strong (Fig. |7|) and weak (Fig. |HJ) coupling data. The Fisher plot of (xx) 2 
versus m Q /(xx) is devised so that data described by a Landau- Ginzburg equation of state 
(ie. 5 = 3 and p = 1 in (fT3*j0 yield trajectories of constant (3 as parallel straight lines, 
intersecting the |/-axis for (3 < (3 C , the x-axis for (3 > f3 C) and the origin for (3 = (3 C . 

Overall, the quality of the fit is remarkably good. The only sign of inconsistency is the 
disagreement in the shape of the curve in the weak coupling phase with (3 > 0.55, although 
even in this regime the fit apparently passes close to "centre of mass" of the data. As seen 
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in Fig. |S1 the data here lie below the fit but curve more strongly, which they must do if they 
are ultimately to intercept the y-axis in the chiral limit. 

For Nf = 4, we were unable to find any fit of comparable quality. Global fits to all 143 
points in the dataset typically yielded a x 2 /dof ~ 0(5). Since in Fig. |2]the low mass points 
are clearly finite volume effected, one strategy is to exclude them and only fit to the 42 
datapoints with (3 < 0.20, m > 0.01, yielding 

A = 2.27(29) 5 = 2.59(18) #. = 0.212(4) X 2 /dof = 0.95 , . 

B = 0.60(7) p = 0.98(13) ^ } 

A slightly more satisfactory fit including a wider window in (3 was found when some attempt 
at incorporating a finite-volume scaling analysis was made, by modifying the equation of 
state to [T7| 

m = A((f3 - (3 C ) + CL^)( XX ) P + B(xx) 5 , (16) 
where for the exponent v we use the hyperscaling prediction 

The best fit to (|16|) was found for the the 121 datapoints with 0.15 < (3 < 0.4, itlq < 0.015: 

A = 1.988(24) 5 = 2.653(30) (3 C = 0.2129(8) x 2 /dof = 8A 

£ = 0.451(9) p= 1.152(3) C = 0.176(7) ( ' 

Although the x 2 /dof is much higher, this fit does a better job in the "symmetric" phase 
(3 > 0.4, and is shown by dashed lines for in Fig. |21 and in Fisher coordinates in Fig. El 
(strong coupling) and Fig. (weak coupling). In any case, the fitted values of the critical 
parameters do not vary much between ()15j) and (|18|). 

As for Nf = 1, the Fisher plot highlights systematic discrepancies in the symmetric 
phase. Overall, though, it is difficult to conclude on the basis of this analysis whether chiral 
symmetry is restored for Nf = 4 as (3 — > oo or not. There is, however, a significant change 
in the fitted exponent f3 mag as Nf changes from 1, where from ()14j) it has the value 1.20(3), 
to Nf = 4 where it takes the value 0.62(9) (fTB^) or 0.67(1) (|18|). This trend is consistent with 
the distinct curvatures of the (xx)o(P) data in Fig. El and is a more dramatic effect than 
the increase of 5(Nf), which even for Nf = 4 still lies below the Landau- Ginzburg value. 



4 Discussion 

Our main result is that there is convincing evidence for spontaneously broken chiral symmetry 
in the continuum limit in the theory with Nf = 1; any extrapolation of the data of Fig. 
resulting in a vanishing condensate in this limit must rely on a hitherto unobserved curvature 
close to the chiral limit. The one theoretical framework in which we might understand 
the origin of such a curvature, namely chiral effective theory, has the assumption of chiral 
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symmetry breaking built in. Therefore we can say with some confidence that Nf c > 1, and 
with slightly less confidence that the dimensionless condensate f5 2 {^ip) ~ O(10~ 3 ). 

For Nf = 4 the situation is much less clear cut. We could not possibly hope to confirm 
Nf c > 4 by direct simulation using current resources, since even for Nf = 2 the dimensionless 
condensate, if there, is too small to measure j^j. We have instead contrasted Nf = 4 with 
Nf = 1 by treating the passage from strong to weak coupling as if there were a true phase 
transition, and attempted to characterise the equation of state by critical exponents. It 
is somewhat ironic that this proved quite successful for Nf = 1 where we believe there is 
no transition but merely a crossover; indeed the quality of the fit is if anything superior 
to similar fits obtained in systems where a true phase transition is believed to occur (eg. 
[T7| |2H1 HE])- For Nf = 4 the fits are less compelling, although this could merely reflect a 
less comprehensive dataset. The only significant difference between the two theories is in the 
exponent f3 mag , which halves as Nf increases from 1 to 4. There is also a small increase in 
S from 2.3 to 2.7, which should be contrasted with the chiral symmetry restoring transition 
in the 2+ld Thirring model (which has the same global symmetries as QED3 in both lattice 
and continuum theories), which for Nf = 2 has 5 ~ 2.7, increasing to 3.1 for Nf = 3, 3.4 for 
N f = 4, and finally to > 4 for N f = 5 [T7HT%1 IT3] (data for the Thirring model with N f = 1 
do not exist). Of course, the hypothesis that the two models have a common fixed point 
behaviour implying the exponents should coincide only holds for Nf = Nf c . Naively one 
might deduce from the trends in these numbers that this would occur for Nf < 2; however, 
one should bear in mind that the equivalence should only hold for Thirring transitions at 
G 2 = 00 and QED transitions at f3 — 00. In the Thirring case there is some ambiguity about 
identifying the strong coupling limit using a lattice regularisation and it is unlikely 
that such information about QED could ever emerge from a numerical simulation, though 
of course, it may be an interesting question to address for analytic approaches. It is also 
interesting to speculate about the nature of the universality class of the transition at finite f3 c 
for Nf > Nf c , in particular the question of whether it is governed by the global symmetries 
of the continuum or the lattice model. 
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Figure 1: Chiral condensate (xx) vs - mass m for Nf = 1 from 24 3 (diamonds), 32 3 (filled 
circles), 44 3 (triangles), 48 3 (squares) and 80 3 (V) lattices. The lines denote fits to data at 
constant j3 to be discussed in Sec. 13. 3( for j3 values (from the top) 0.25, 0.275, 0.3, 0.325, 
0.35, 0.375, 0.4, 0.45, 0.55, 0.65 and 0.90. 
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Figure 2: (xx) vs - m o f° r -W/ = 4 from 16 3 (triangles), 24 3 (diamonds), 32 3 (filled circles) 
and 48 3 (squares) lattices. The lines denote finite volume scaling fits to data at constant (3 
(see Sec. USD, for (3 values (from the top) 0.15, 0.1625, 0.175, 0.1875, 0.2, 0.2125, 0.225, 0.4 
and 0.6, with dotted lines denoting fits to 16 3 and dashed lines fits to all other volumes. 
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Figure 3: The value of (xx) obtained in the chiral limit vs. (3 
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Figure 4: /3 2 (xx) vs - /3 m o for iV/ = 1,2 and 4 for the datasets closest to continuum and 
chiral limits. 
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Figure 5: (3 {xx) vs - P m o fo r -W/ = 1 in the region of the chiral limit. 
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Figure 6: Dimensionless masses (3M n vs. (3 for the various pion operators discussed in 
Sec. EH 
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Figure 7: Fisher plot (xx) 2 vs - m o/(xx) f° r strong coupling data with Nf — 1. Unless 
stated all data are from 24 3 lattices. The dashed lines are the result of the fit (|14|). 
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Figure 8: Fisher plot for weak coupling data with Nf = 1. 
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Figure 9: Fisher plot of the strong coupling data from a 24 3 lattice for Nf = 4; the dashed 
lines are the result of the fit (jTSjl . 
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Figure 10: Fisher plot of the weak coupling data for Nf = 4. 
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